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Abstract 

We consider the gravitational generation of the massive Z boson field of the standard model, due to 
the natural breaking of its conformal invariance during inflation. The electroweak symmetry restoration 
at the end of inflation turns the almost scale-invariant superhorizon Z spectrum into a hypermagnetic 
field, which transforms into a regular magnetic field at the electroweak phase transition. The mechanism 
is generic and is shown to generate a superhorizon spectrum of the form B a 1/f on a length-scale £ 
regardless of the choice of inflationary model. Scaled to the epoch of galaxy formation such a field suffices 
to trigger the galactic dynamo and explain the observed galactic magnetic fields in the case of a spatially 
flat, dark energy dominated Universe with GUT-scale inflation. The possibility of further amplification of 
the generated field by preheating is also investigated. To this end we study a model of Supersymmetric 
Hybrid Inflation with a Flipped SU(5) grand unified symmetry group. 
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1 Introduction 



Magnetic fields permeate most astrophysical systems and their presence may have numerous 
cosmological and astrophysical implications. Indeed, magnetic fields may substantially influence 
the formation process of large-scale structure || and of individual galaxies However, the 

most important role of large-scale magnetic fields is that they may be responsible for the magnetic 
fields in galaxies. 

It is a well-known observational fact that galaxies feature magnetic fields of strength ~ /iGauss 
Hi H • The structure of such fields in spiral galaxies follows closely the spiral pattern |5| and 
this strongly suggests that the galactic magnetic fields are generated and sustained by a dynamo 
mechanism According to the galactic dynamo, the cyclonic turbulent motion of ionized 

gas combined with the differential galactic rotation amplifies a weak seed field exponentially until 
the backreaction of the plasma motion counteracts the growth of the field and stabilizes it to 
dynamical equipartition strength. However, the origin of the required seed field remains elusive. 

In order to trigger successfully the galactic dynamo, the seed field has to satisfy certain re- 
quirements of strength and coherence. Indeed, it has been shown that seed fields which are too 
incoherent may destabilize and destroy the dynamo action J8| . Most dynamos require a minimum 
coherence length comparable to the dimensions of the largest turbulent eddy, ~ 100 pc. The 
required strength is determined by the dynamo's amplification timescale (typically the galactic 
rotation period) and the age of the galaxy. Recent observations suggest that the Universe at 
present is dominated by a dark-energy component || . In such a case the galaxies are older than 
previously thought, and the minimum seed-field strength may be as low as -B S ccd ~ 10 _30 Gauss 

There have been attempts to generate the necessary seed field using various astrophysical 
mechanisms, the most important of which involve battery |llj or vorticity effects Battery 
mechanisms require a large-scale misalignment of density and pressure (temperature) gradients, 
usually associated with large lobe-jets (AGNs) or starburst activity Jig ], and are therefore difficult 
to realize in the majority of the galaxies. On the other hand, large-scale vortical motions can be 
effective only if the ionization of the plasma is substantial, which can hardly occur as late as the 
epoch of galaxy formation. 

Due to the above difficulties, it has been frequently argued that the origin of the seed field may 
precede the galaxies themselves, i.e. be truly primordial. The obvious requirement for large-scale 
magnetic-field generation before the time of recombination is that it has to occur out of thermal 



equilibrium, because such a field breaks isotropy |14 . This limits the choice to generating the 
magnetic field either at a phase transition or during inflation, although more exotic mechanisms 
have occasionally been suggested, such as vorticity-inducing cosmic strings or magnetic fields 
in the pre-Big Bang era Jig] . 

There have been numerous attempts to create a primordial magnetic field at the breaking of 
grand unification or at the electroweak phase transition or even at the quark-confinement epoch 
(QCD transition) 

However, since the generating mechanisms are causal, the coherence of 
the created magnetic field cannot be larger than the particle horizon at the time of the phase 
transition. Because all the above transitions occurred very early in the Universe's history, the 
comoving size of the horizon is rather small (the best case is the QCD transition, for which the 
horizon corresponds to ~ 1 AU) and so the resulting magnetic fields are too incoherent. It has 
been shown that not even the most favorable turbulent evolution can adequately increase the 
correlation length of such fields jl8|. Thus, one should achieve super horizon correlations in order 
to generate a sufficiently coherent magnetic field in the early Universe. 

Inflation presents the only known way to achieve correlations beyond the horizon scale and, for 
this reason, has received a lot of attention. However, the prime obstacle to generating magnetic 
fields during inflation is the conformal invariance of electromagnetism, which forces the magnetic 
flux to be conserved fljfl . As a result, the strength of any generated magnetic field decreases 
exponentially as the inflationary Universe rapidly expands. Attempts to overcome this problem 
have set out to break the conformal invariance in various ways, such as through the explicit 
introduction of terms in the Lagrangian which couple the photon directly to gravity or to a 
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scalar or pseudoscalar field, or through the inclusion of a massive photon, or even by means of 
the conformal anomaly p^|[pO[. However, those very few attempts that succeed in producing a 
sufficiently strong magnetic field do so relinquishing simplicity. Recently, interesting efforts have 
been made to create magnetic fields by coupling the photon to some scalar field during inflation 
pf[ or at the preheating stage 22 1 . These proposals have since been criticized in J23| . Other recent 
ideas include magnetogenesis due to the breakdown of Lorentz invariance in the context of string 
theory and non-commutative VSL theories J24| , due to the dynamics of large extra dimensions pq ] 
and, finally, due to gauge field coupling to metric perturbations p6| . The viability of the latter 
has been questioned in |27|] . 

Recently one of us (TP) has studied the effects of conformal symmetry breaking due to the 
coupling of the photon field to fermions and scalars [^7|. The author has considered effective 
actions arising from loop corrections in the 1/M expansion and from the anomaly. In addition, 
the photon coupling to scalar and pseudoscalar fields was reconsidered, having been originally 
introduced by Ratra in (20) and in |l^] (see also Garretson et al. in |^0| and |2^|). 

In this paper we show that conformal invariance is naturally broken in inflation without the 
need for any exotic mechanism or field. In particular, we consider the gravitational generation of 
the massive Z-boson field of the standard model. The electroweak symmetry restoration at the end 
of inflation turns the superhorizon Z spectrum into a hypermagnetic field, which transforms into a 
regular magnetic field at the electroweak phase transition. The mechanism is rather generic and is 
shown to generate a superhorizon spectrum of the form B lms cx l/£ on a length-scale t regardless 
of the choice of inflationary model. However, it is possible to consider further amplification of 
the generated field via parametric resonance during preheating. To this end, we study a model of 
Supersymmetric Hybrid Inflation (SUSY-HI) with Flipped SU(5) as the symmetry group of the 
Grand Unified Theory (GUT). Our mechanism for creating the field has been described in [ p9[ , 
where it was not specifically applied to the Z field. 

The structure of our paper is as follows. In Section 2, the model of Supersymmetric Hybrid 
Inflation is presented. Section 3 contains a study of the relevant particle representations of the 
Flipped SU(5) group. The symmetry breaking process is analyzed and the field equations of 
the gauge bosons are obtained, with particular emphasis on the hypercharge field and its source 
current. Also, the equations for the scalar fields (Higgs fields and the inflaton) are layed out in 
detail. In Section 4, we study the gravitational production of the Z-boson field during inflation in a 
model-independent way. The resulting spectrum of the magnetic field is computed and scaled down 
to the epoch of galaxy formation, thus showing that it may be sufficient for explaining the galactic 
magnetic fields. In Section 5, the possibility of extra amplification by preheating is investigated 
using the particular model of Flipped-SU(5) SUSY-HI. Both Flipped SU(5) and SUSY-HI are well 
motivated by Supergravity and Superstrings and rather are naturally compatible. However, it 
should be noted that the choice of GUT group was made only to allow us to perform analytical 
calculations. In principle, any GUT group would suffice. Finally, in Section 6 we present our 
conclusions. We have attached two Appendices, one describing the Hartree approximation used at 
certain points in the calculations and the other presenting the way to calculate the rms amplitudes 
of the Z-boson and the hypermagnetic field. 

Throughout the paper we use a (+, — , — , — ) metric and units with c = h = 1 so that Newton's 
gravitational constant is 8nG = m^p 2 , where nip = 2.4 x I0 18 GeV is the reduced Planck mass. 



2 Supersymmetric Hybrid Inflation 
2.1 Hybrid Inflation 

Hybrid Inflation was originally suggested by Linde |]3(][| to avoid the fine-tuning problems of most 
inflationary models due to radiative corrections. In order to achieve this, Hybrid Inflation intro- 
duces a mass scale M, much smaller than the Planck mass, that sets the scale of the false- vacuum 
energy towards the end of inflation. Thus, at least near the end, the inflationary potential is 
protected from radiative corrections and is sufficiently flat. The expense paid is the necessity of 
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introducing, in addition to the inflaton field s, another scalar field (ft related to the scale M. The 
scalar potential for Hybrid Inflation is, 

V((ft, s) = -A {(ft 2 - M 2 ) 2 + l -hs 2 4> 2 + V 3 (s) (1) 

where A and h are coupling constants and V s (s) is the slow-roll potential for the inflaton field. 
Because of the coupling between s and (ft, for s > s c , where s c = M^J X/h, the above potential 
has a global minimum at (ft = 0, so that V(0, s) = AAf 4 /4 and we obtain low energy-scale inflation 
as required. However, when s < s c , spontaneous symmetry breaking displaces the minimum of 
the potential (||) to (ft — M e ff, where M^t = M 2 [l — (s/s c ) 2 ]. The system rapidly rolls towards 
the new minimum and oscillates around it, thus terminating inflation. In general, inflation ends 
abruptly less than one e-folding after the symmetry breaking. 

Letting M be of the order of the grand-unification scale can indeed satisfy the requirements 
of successful inflation, for which the total number of e-foldings must be large enough (typically 
about 60) to solve the flatness and horizon problems and to account for the magnitude of density 
perturbations and the anisotropy of the Cosmic Microwave Background Radiation (CMBR) in 
agreement with Large-Scale Structure and COBE observations. Thus, a natural candidate for the 
(ft scalar is the Higgs field of a Grand Unified Theory (GUT). Such a choice has the merit of not 
introducing any additional unknown scalars or mass scales. 

However, the most important attribute of Hybrid Inflation is that it can originate from Super- 
symmetry and is, therefore, one of the few inflationary models with a theoretical foundation in 
particle physics. 



2.2 Supersymmetric model 

The literature on Supersymmetric Hybrid Inflation (SUSY-HI) is rather rich, as it is possible to 
attain from either Supergravity or Superstrings pi] . We shall briefly describe an F-term GUT 
inflationary model as in p2| . D-tcrm models also exist p3| , but will not concern us here. 
The most general renormalizable superpotential with i?-symmetry is p2| , 

W = kS(<P$ - M 2 ) (2) 

where "J"! is a conjugate pair of singlet components of chiral superfields that belong to a nontrivial 
representation of the GUT group G, S is a gauge-singlet superfield, and k, M are constants that 
can be made positive by phase redefinitions. Introducing the above expression to the _F-tcrm 
scalar potential Vf — Y^i \dW/d(fti\ 2 one finds, 

V = k 2 \M 2 - $$| 2 + n 2 \S\ 2 (m 2 + |$| 2 ) (3) 

where we write the scalar components with the same symbols as the superfields. Using i?-symmetry 
we can bring the scalar fields onto the real axis. We set S = and $ = $ = (ft, where s and (ft 

are real scalar fields]]] Then the scalar potential becomes, 

V((ft,s) = n 2 {(ft 2 -M 2 ) 2 + k 2 s 2 ^ (4) 

This is, in fact, the potential for Hybrid Inflation given in (0) with \ = 2h= 4k 2 and s c = \[2M. 

The non-vanishing vacuum energy breaks supersymmetry and generates radiative corrections, 
which induce the required slow- roll potential V s (s). It can be shown that the overall contribution 
is of the form |32| , 

where A is a suitable renormalization mass scale. Thus, the radiative corrections provide a gentle 
down-slope for the inflaton, V s oc In s, which helps to drive it towards its minimum. 



The D-term vanishes because <E> = <3? corresponds to the D-flat direction that contains the supersymmetric 
vacua. 
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At earlier stages of the inflaton's evolution, s may be of the order of the Planck mass, so that 
Supergravity corrections should also be considered. However, it can be shown |34j that the flatness 
of the potential is preserved for minimal Kahler potential (K = 2(f) 2 + S 2 ). 



3 Flipped SU(5) 

In choosing the GUT for the symmetry breaking that terminates Hybrid Inflation, one has to 
ascertain that there is no monopole problem |]35| . Thus, a simple or semi-simple group is not 
an option. We decided in favor of Flipped SU(5) because of its simplicity and its resemblance 
to the Standard Model (SM). Flipped SU(5) is jargon for the GUT group SU(5)xU(l), in which 
the hypercharge U(l)y is contained both in the SU(5) part and the U(l), in contrast to the 
Georgi-Glashow SU(5)xU(l) model, which has the hypercharge fully embedded in SU(5) [36| . 

The supersymmetric version of Flipped SU(5) is well motivated by Superstrings Moreover, 
it can be considered as an intermediate stage in the breaking of supersymmetric SO(10) J3(| |38| : 



SO(10) -> SU(5)xU(l) -> SU(3) c xSU(2)xU(l) y 

Since the breaking of SO (10) to Flipped SU(5) generates monopoles, it would have to take 
place before or during the inflationary period, so that the monopoles can be safely inflated away. 
The SUSY-HI in this model has been studied in [||. 

3.1 The model 

The Lagrangian density of the model is , 

£ = -\g^g^MG^G p<T ) - \g w g v °G%G* pa + l -g^Tr[{D^D^] + ^ (V^V.s - V (6) 

where $ is the GUT-Higgs field, s is the inflaton field, G pil and G° are the field strengths of the 
SU(5) and the U(l) gauge fields respectively, is the covariant derivative of the Higgs field, 

and V is the effective potential. 
The SU(5) field strength is, 

G M „ = V P A V - V V A„ - iglA^Au] (7) 

where g is the gauge coupling of SU(5), A p = A^T a , and A% (a = 1...24) are the SU(5) gauge 
fields with T a being the corresponding generators. Similarly, for we have, 

G% = V^Al - V v Al (8) 

where A^ is the Abelian gauge field of the U(l) group. 

Flipped SU(5) is broken by the GUT Higgs field $, which belongs to the (10,1) antisymmetric 
representation. The U(l) degree of freedom corresponds to an overall phase. In this representation, 
the generators T a are given by T a — ~^fi a , where the set of modified 5x5 Gell-Mann matrices 

fi a can be found in the appendix of jh|. Also, for the Abelian generator we can write T° = 

and /i° = *J\l with / being the 5x5 identity matrix. The normalization constant has been chosen 
so as to simplify the treatment of the symmetry breaking process. 

Writing G^ = G a ^T a and A p = A a ^T a , the field strength (g) becomes, 

G% - d^Al - d v Al + gr hc A\Al (9) 

where we have used that [T a ,T b ]= if abc T c and Tr(T a T b ) = S ab with f abc being the structure 
constants of the group. Thus, in (g) we have Tr(G ^G pa ) — G a ^ v G a pa . 
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The covariant derivative is, 



= - igd{A^ $) - igd(A^ $) 
where g is the U(l) gauge coupling, A^ — yl°T , and, 

Since <£>,j is antisymmetric, the above can be written in matrix notation as, 

d{A^ $) = A^ - (A^f = ^[T a $ - (T Q $) T ] 

and, similarly, d(A fi ; $) = -^6/5 71°$. 

The matrix of the gauge fields, which is hermitian, may be written as, 

-V2(A U . + ^A^ 



(10) 



(11) 



(12) 



n 
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/G 2 H 




c'yV + CyY 


G 1 




G 3 


X 1 


F 1 \ 




G 1 


-G 2 + 




Yc Y Y 


G 4 


x 2 


F 2 




G 3 




G 4 






-c Y Y X 3 


F 3 




X 1 




X 2 






w 3 - 




V 


y 1 








F 3 







(13) 



where we have suppressed the spacetime indices inside the matrix. In the above we have employed 



the following definitions: Q 










XI = Al + 


iAf 




G\ = A\ + 


iA\ 


*2 = 4 1 - 


M< 




^ = 4 




XI = Af - 


H< 




G 3 = 4 + 


iA\ 


Y^=Afl 


-< 




g^ = 4 + 


iA 7 , 


Yl = A? H 






^ = 4 












W~ = yl 22 


+ ^ 3 




w* = 4 4 




where 












-VI 


fl 


i 

Cy = Cy 

Cy 


9v 


and also, 










2 In what follows in this Section the 


over-bar denotes 


charge conjugation [e.g 





(14) 



9 9 

g 2 + g 2 



(15) 



U(l). 
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9 A ° 
9 



(16) 



In the above, the complex bosons X p , Y? (a = 1, 2, 3) are the supermassive GUT-bosons, 
are the massless SU(3) C bosons, W p and W p are the usual SM VU-bosons and Y p is the hypercharge 
gauge boson, which should not be confused with the massive Y^ s. Finally, is a massive boson 
generated partly by SU(5) and partly by U(l). It can be viewed as the analogue of the Z-boson 
of the SM. 

In analogy with the matrix, one can define the matrix 



H, 



lv — — V2 [^g^v 



^ fJ,L> 



Using this and (0) it is easy to show that 



(17) 



(18) 



Without loss of generality the GUT Higgs field may be chosen to lie in the following direction: 



V2 



( 











\ 












































1 


V o 








-1 





(19) 



where <p is a real positive function of time <f> = 4>(t). In this case, it can be shown that the 
interaction term of the Lagrangian density becomes, 



£ ln t = -g^Tv[(D^yD^} = -g^{d^){d v <f>) + -M' XY g^ {XpT v + Y^) + ~-M v g^V v 

(20) 

(21) 



where the masses of the supermassive GUT-bosons and of V p are 



M XY = -i 



V2- 



M v = — M XY 
c v 



Thus, the interaction term reduces to the mass terms of the massive bosons plus the kinetic 
term of cf>. Note that after the GUT symmetry breaking the W's and the hypercharge K, remain 
massless. 

It is straightforward to show that the kinetic term of the supermassive GUT bosons is, 



nXY 
•'-'kin 



4 J J —pa , 

where0iJ^ a = {H^) 4a and = {H^) 5a with a = 1,2,3 



(22) 



3.2 The field equations of the gauge bosons 

The field equations of the massive GUT-bosons are, 



[d M + (d h In y/=D g )] g^g v "{d p X: - d a X«) + M XY g^X« = J, 



[d„ + (d M In y/=D g )] g^g^(d p Y a a ~ d„Y p a ) 



M 2 XY g^ = 



^Because the matrix is also hermitian, we have: Tr(H^ u Hp a ) = {H^ u )ij{H p(7 )i 



5<v 



(23) 
(24) 



G 



where D g = det(g^) and the current, in matrix notation, is,0 



The field equations for the massless VF-bosons are, 

[0„ + (0 M In y/=D^] \g^g^(d p W+ - d CT ^ p +) 



45 



and 



1 



+ (0„ In v^)] - 9 CT % 3 )J = 2^44 ~ J&) 

The field equation for W~ is just the complex conjugate of (p6|). 
The field equations for Y" M and are, 



(25) 



(26) 



(27) 



[fy + (fy ^ v^)] - dvVp)] + M v g^V p 



— COS Jy 



(28) 



\d p + (s M in y^n^v^p^ - W] 



(29) 



where tan 9 = g/g is the GUT equivalent of the Weinberg angle and the hypercharge source 
current is, 



l v — l v 4- l v 



(30) 



3.3 More on the hypercharge source current 

It can be shown that the hypercharge source current can be written as, 

J Y = V2 ff Im{[0 M + (d M m (3i 



The above gives, 



J Y = V2glm{ [8 P + (d p In ^D g )\g^ g v ° {X a p X a a + 
+ g^g^[X^d P X + Y p a d p Y a a - (p ~ a)]} + 



(32) 



+ ffW T Re[(^)4fe(np)fe;(o ff ) i 4 + (nMn p )ki(n a )i5 - (p «-» a-)] 



It is now evident that, only through the last term in the above expression, may the hypercharge 
source current receive contributions from gauge fields other than the supermassive GUT-bosons. 
It can be shown that the contributions to J Y from the various gauge fields are of the following 
form, 

The contribution from the hypercharge field: 



J Y [Y] = g 2 c Y gWg™[Y p Re(X«X a + Y«Y a ) - (p «-> a) 
The contribution from massive V: 

J Y [V] = g 2 c V 9 w 9 ua \ V p Re(X*X a a + Y«Y°) - (p «-> a) 



(33) 



(34) 



4 Note that the current matrix is also hermitian, JY. = Jj t . 
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The contribution from W 3 : 



J Y [W 3 



W 3 Re(Y«Y a a - X«X«) - (p ^ a) 



The contribution from W : 

JyW*] = -g 2 g^g^Re 
The contribution from the gluons G l : 



Y?X a a W+ + X°Y°W p -(p~a) 



where a, (3 = 1, 2, 3. For a — (3 one should consider only the G-part of O"^. 
In view of all the above we can write the hypercharge field equation as, 



(35) 



(36) 



(37) 



[d, L + (0„ In ^D- g )W p 9 ua (d P Y« - d a Y„)\ + f^^Rep^ + *p 
= - 1 cot 6 g 2 Y £^Re(X«X a a + Y«Y*) V p - 

'g Y lm{[d p + (d p In J=E?j\g>>'> (XfX* a + Y«Y^) + 

+ gWg™[X«d p xZ + Y°d p Y a a - (p «-> a)}} - 



g Y g£^P° 



W 3 Re{Y«Y" a -Xpr„)-Re{Y«XZW++X«Y';W-) -J± sin9J£[G] 



(38) 



where we have used that, 



gVPg-** 



A a B p -(p^a) =S^A a B p £^ pa = e reA ^VV 4 



(39) 



for any A a and B p , where e K/iI/ is the totally antisymmetric Levi-Civita tensor. 

Although the above equation appears rather complicated it can be understood as follows. If 
we ignore the expansion of the Universe then (p8[) becomes, schematically, 



□Y + g Y (XX)Y = - cot 9 g Y (XX)V - g Y {(XdX) + (g/V2)[(XXW) + (JIG)]} 



V = Ia [Tr(*f $) - M 2 ] 2 + i^s 2 Tr($t$) + v s (s) 



(40) 



where we have symbolized all the supermassive GUT-bosons with X, we have setj^jgy = a/5/6 g Y 
and we have taken the equivalent of the Lorentz gauge for the hypercharge held. 

3.4 The field equations for the scalar fields 

The effective potential is taken to be, 



(41) 



where M is the scale of the GUT symmetry breaking, A and h are coupling constants and V s (s) 
is the slow-roll potential for the inflaton field. For the particular <I>-gauge introduced in ( |l9| ) the 
scalar potential reduces to the one given in ([!]), for which in SUSY-HI we have A = 2h ~ 1 and 
V s oc In s. 

The field equation of the inflaton is, 



[d p + (0 p In ^D~)] gTduS = ~ 



(42) 



1 Corresponding to the definition used in ]36|. 
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For a spatially- flat Friedman-Robertson- Walker (FRW) metric and with s = s(t) the above 
reduces to the well known form, 



dV 

a + 3Hs + — = 

OS 



(43) 



where H = a/a is the Hubble parameter with a = a(t) being the scale factor of the Universe and 
the dot denotes derivative with respect to the cosmic time t. For the potential (Q) we have, 



dv 2 av. 

— =h(j)S+ — 

os os 

Using a Conformal-time FRW (CFRW) metric and with V s oc In s we find, 



— a + [hip a + a —— j 
a oa 



a" a + (hip"<T + a'— — ) = 

where V a oc lncr, the primes denote derivatives with respect to conformal time t and, 



(44) 



(45) 



a — a s if = a<j) (46) 

The field equation for $ is more complicated. The general expression is found to be, 



[0„ + (d„ In y/=D~ g )\ gT&v*) = igg^A;d(T a ; £>„$) - {a [Tr(&§) - M 2 ] + hs 2 }$ (47) 
The above is a complex matrix equation. The resulting relevant equation for (p in CFRW is, 



1 



„w« 4 



^ -—<p+ [A(^ - a 2 M 2 ) + ha 2 } <p = -a 2 g 2 g^(X«X„ + Y«Y v + —V^V V ) <p 
where we ignore all but the spatially uniform mode of ip. 



(48) 



3.5 The contribution of the electroweak Higgs field 

Due to supercooling (T a a -1 — ► 0) the electroweak (EW) symmetry is broken during inflation 
and, thus, it is important to determine the contribution of the EW-Higgs field W to the field 
equations of the gauge fields. To that end we have to consider the embedding pattern of the SM 
group into our GUT. 

The SM group is not fully contained in the SU(5) part of Flipped SU(5). The SU(2) and the 
SU(3) C of the SM are contained in SU(5) so that their couplings, say g<i and 33 merge and equal 
g at the GUT scale. But the hypercharge coupling gy does not. In fact this coupling can be 
considered to be comprised of the g coupling of U(l) and a g\ coupling corresponding to a U(l)i 
subgroup of the SU(5) |3(|. This 171 merges with 172 and g-$ at the GUT scale. Thus the structure 
of the symmetry breaking is, 



Flippcd-SU(5) 
SU(5) -> 



SU(3)xSU(2) 



U(l) 
where 




U(l)y 



Standard Model 



SU(3)xSU(2)xU(l)i 



91 

.92 

93 



at GUT scale 



and 



gr = 



fjfj 



VgJW 2 
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It is obvious that at the GUT scale, because g\ — > g, the hypercharge coupling is given by ( |l5| ) 
as required. The U(l)i generator corresponds to the T 15 generator of SU(5) since this is the one 
that gets mixed with the generator of U(X) to give the hypercharge generator. 

In our framework we have the peculiar situation that, although the GUT symmetry is un- 
broken during inflation due to the coupling between the inflaton and the GUT-Higgs field, the 
EW symmetry is, in principle, broken because the Universe is supercooled.^ Therefore, we are 
interested in finding the contributions of the EW-Higgs are to the field equations. 

The additional EW contribution to the Lagrangian (||) is of the form, 

Cbw = \g^{D^D v ^ - Cf(*t*) (49) 

where, 

f/(*t$) = ^K{^^ - M 2 EW f (50) 



with Mew being the electroweak scale and 

= (V„ + ig2W«r a + igyYp)* (51) 

where r Q are the generators of the SU(2) of the SM such that Wj}r a = -j2A° il T a with a = 21 + a, 
a = 1,2,3 and at GUT scale g2 — > g- 

The EW-Higgs field is in the 5 complex vector representation. Because the gauge condition 
( |l9| ) leaves the full EW symmetry group unbroken, we can, without loss of generality, introduce 
the gauge, ^f T = ip(0 1 0), where ip — ip(t) is a real, positive function. Then the complete field 
equations of the Ws are found to be, 



[d, + (d. In y/^)] [gWg™(d p W+ - d a W+)\ + M w g^W+ = J£ 5 (52) 

and 



[d p + {d^\n^jT g )] g» p g v °{d p Wl-d a Wl) +M w g^Wl = -(jr 4 -J 5 ^)+sine( 5 V)V% (53) 



where Mw = gij) and the field equation of W„ is just the complex conjugate of (]5 
Moreover for the hypercharge we have, 



ld p + (d p \n^D~ g )}{g^g™(d p Y a ~d a Y p )} + M 2 g^ = - /J sin 9 J y +sin9(<^) VWjJ (54) 

where My = gyip- It can be shown that the contribution of the EW-Higgs to the field equation 
of the massive boson is zero, as expected. Now, let us rotate the above to form the Z p boson 
and the photon A p of the SM. They are defined as, 



= cos6 w Wl-su\0 w Y p (55) 



A p = sm6 w W p + cos 6 w Y p (56) 



where 9w is the Weinberg angle defined as, tan 8w = gy /<?2 • Note that at GUT scale tan 6w — sin 9 . 
Thus, we obtain, 

6 In fact, the quantum fluctuations restore also the electroweak symmetry in a se nse. What actually happens is 



that the EW-Higgs field forms a non-zero condensate, as we will explain in Sec. 4.3, which provides masses for the 
massive EW gauge bosons, thereby breaking their conformal invariance. 
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[0 M + {d, In y/^DHyVfridpZ. - d a Z p )\ + Mf 5 ^ 



2 cos t% 



(1 + sin 2 0)J& - (1 - J§ sin 2 6)^ 



where 



and 



COS t/^y 

Thus, we see that Mz = Myy/ cos 8\y as expected. For the photon we find, 



(57) 



(58) 



[d„ + {d^^D^Wg^idpAv - d a A p )] = isin%[(l 
Finally, the field equation of ^ is, 



|)J44-(1 



3 / J S5 



(59) 



[fy + (fy In v^)] (/""(A,*) = _^( 5W -a T a + 9r Y^)£>^ _ A,(*t$ _ m 2 ew )V (60) 
which, in the gauge mentioned above, becomes, 



where 



y" --y+ [K(y 2 a 2 M% w )]y = a 2 g^[g\WlWl + W 2 W 2 U ) + g 2 z ZpZ„} y (61) 



y — aip (62) 



It should be pointed out here that the gauge boson masses appearing in the corresponding field 
equations above are a purely mathematical result of the structure of the Flipped SU(5) group. 
They will be non-zero only if a Higgs-field condensate is generated, such as happens in the relevant 
phase transition. 



4 Conformal invariance breakdown during inflation 

In this section we will show that the Z-field is naturally produced during inflation because its 
mass term is non-zero and, therefore, the field is not conformally invariant, but instead it has 
gravitational source terms. At the end of inflation reheating restores the EW-symmetry and 
the generated Z-spectrum is projected onto the direction of the massless, Abelian hypercharge 
field, thus, creating a hypermagnetic field, which freezes into the primordial plasma and survives 
until the EW phase transition. At that time the hypercharge configuration is projected onto the 
photon, transforming the hypermagnetic field into a regular magnetic field, which evolves until 
galaxy formation. We will show that such a field may be strong enough to seed the dynamo in 
galaxies and account for their observed magnetic fields. 

4.1 The mode equations for Z 

During inflation the temperature is essentially zero, which means that the electroweak symmetry 
is broken. Thus, we expect tp ^ 0, which, in view of ( |58| ) renders the Z^ gauge field massive. On 
the other hand, due to the interaction between the inflaton and the GUT-Higgs field, the GUT 
symmetry is unbroken and <f> — 0. This means that the GUT-bosons X" and Y" are massless. 
The same is true also for the bosons because they do not couple to any scalar field. Therefore, 
the X^, Yjf and the bosons remain conformally invariant during inflation so, like the photon, 
they cannot be generated gravitationally. Thus, their magnitude during inflation is negligible. For 
this reason, because the current in (|57|) is primarily sourced by the GUT-bosons we will ignore its 
contribution during inflation. 
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Assuming a CFRW metric <j^„ = c^r) 2 ^^ the field equation ( p7[ ) may be rewritten as, 

rff>rf" , d ll {d p Z a - d a Z p ) + a 2 M 2 V ^Z p = (63) 

where rf v is the flat spacetime Minkowski metric and Mz = Mz(t) is the mass of the Z-boson 
generated primarily by the self-interaction term of the EW Higgs field during inflation as will be 
discussed in more detail in Sec. 4.3. 

The temporal and spatial components of the above give, 



<9 T (V • Z) - V 2 Z T + a 2 M 2 z Z T = (64) 

UZ- V(d T Z T - V ■ Z)+a 2 M 2 Z = (65) 

where V 2 = didi with i = 1, 2, 3 is the Laplacian, □ = d 2 — V 2 is the D'Alembertian and V is the 
divergence or the gradient. Taking the derivative of (p3) we obtain the integrability condition, 

d T Z T - V ■ Z = -2 [d T \n{aM z )]Z T (66) 
in view of which we can recast ( |65| ) as, 

[□ + (aM z ) 2 ]Z + 2 [d T ln(aM z )]VZ T = (67) 
We now switch to momentum space by defining, 

d 3 k 



Z i1 {x,t)= J ^-^Z A1 (fc,r)exp(ife • x) 
Then ( |64|) and (|62]) transform into, 



(68) 



id T (k ■ Z) 

k 2 + (aM z ) 2 = ° (69) 

[d 2 T + k 2 + (aM z ) 2 }Z + 2i[d T \n(aM z )}kZ T = (70) 
which, when combined, result in, 

[d 2 + k 2 + [aM z f]Z + 2 -^f^0 d T [k(k -Z)]=0 (71) 

where k 2 — k • k. We can decompose the gauge field into longitudinal and transverse modes in 
the manner, 

zll ^ k(k^Z) Z ± ^Z- ZW (72) 

Then (|l|) gives, 

{^ + 2 ^+sr » +ti+Mjz)a } z " - ° (73 » 

r d 2 + k 2 + (aM z ) 2 ]z 1 - = (74) 

In the following we will concentrate on the transverse component since the longitudinal com- 
ponent is really relevant only when interactions are taken into account, which, however, are not 
dealt with in this section. For simplicity, we will drop the symbol. 
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4.2 Gravitational production of Z bosons 

In CFRW coordinates the scale factor during inflation and radiation domination is given by, 

a(r) = -1/Ht -oo<t<-1/H 

(75) 

o(t) = Ht t>1/H 

such that a(-l/H) = a(l/H) = 1 and a'(-l/H) = a'(l/H) = H, where H ~ const, is the Hub- 
ble parameter during inflation. Note that we have assumed a sudden transition from inflation 
to radiation. If there is an intermediate era evolving like a matter-dominated Universe this can 
be taken into account; but we do not expect that the conclusions reached below will be affected 
by that, in that the spectrum should remain unaltered. Furthermore, because we are consider- 
ing a GUT-scale inflationary model, we expect that reheating restores the EW symmetry and, 
therefore, Mz(t> H~ v ) = 0. In view of the above the equation for the transverse component of 
Z k (r) ee \Z J -(k,r)\ given by @ becomes^ 



M 2 1 \ 

fc 2 + -^f^JZ fc (r) = (inflation) (76) 
(d 2 + k 2 ) Z fe (T) = (radiation), (77) 



As we will show in Sec. 4.3 , M§ is created by the self- interaction of the EW Higgs-field. Because 
of this there is a slight logarithmic correction due to the time dependence^] of Mz- However, we 
expect that this change does not affect in any significant manner the results presented below. 

Equations (|76| ) and (|77]) have solutions in terms of Bessel and harmonic functions, respectively. 
In particular ( |76| ) may be written as, 

d 2 T + k 2 - V —^jZ k {r)=Q (78) 
with 



which has solutions in terms of the Hankel functions H u , which have the Wronskian normalization 



1 Ml 

4 - W ^ 

rOO 

-v I 

W[Hi 1 \z),Hl 2) (z)} = HP(z)d z Hi 2) (z) - Hl 2) {z)d z H^ ) (z) = -Ai/[irz). Hence, as a basis for the 
vector-field mode functions during inflation, we shall choose, 

ZP(kr) = ly/=mVHP(-kT) (j = 1, 2), (80) 
which solve ( [76]) and have the Wronskian normalization Ell 

W[Zi l \kr),Zl 2 \kr)] ee Z^(kT)d T Z^\kT) - Z^{kT)d r Z^{kr) = Vi (81) 
Asymptotic forms of these solutions are easily found Ell : 



Zy\kr^-oo) ~ y/V/2k e -*(-iy[kT+%(v+i)} v __ + 0[{-kr)- 2 ] \ {j = 1,2) 

(82) 



and also 



7 Since the coefficients in ( |74| ) for Z 1 - do not depend on direction, it is reasonable to assume that (on average) 
the dependence is only on the magnitude of the momentum, k = \k\. 
8 In fact, 1/t 2 becomes — ln(— Ht)/t 2 because M| oc — ln(— tH). 
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>(T) ~ -V-tttV 
; 2 



i(-iy 

2-1/1 



2 

fer 



1 - (-l)^COt(W) 



rr \ /cry r(f + 1) 

Re[v] > (j = l,2) 



Y 



+ 0[(-fcr) 

In the radiation era the eigenfunctions that solve ( |77| ) are 

^(fcr) = N /V/2fc S-V* (i = l,2) 



(83) 



(84) 



which are appropriately normalized as in (|3l|). 

It is now obvious that there is particle production since the eigenfunctions in inflation and 
radiation era are not orthogonal. A simple way of computing particle production is to do the 
matching at the end of inflation. Choosing for example Z^F 1 in inflation, we can match it to the 
radiation era eigenfunctions as follows 



ZP(-k/H) = a k zll ] d (k/H) + (3 k Z^ d (k/H) 
d T Z^ (-k/H) = a k d T Z^ (k/H) + (3 k d T zlH (k/H) 



(85) 



where a k and (3 k are the Bogoliubov coefficients. The asymptotic form of Z^ in ( |S2|) implies that 
the Bogoliubov transformation is diagonal in momentum space, and that the coefficients depend 
only on the magnitude of the momentum k = \k\. A similar transformation can be constructed 

(2) 

for Zu ■ The fact that f3 k ^ implies that there is particle production in inflation. The number 
of particles produced per mode is then of the order of n k ~ \{3 k \ 2 - Using the asymptotic form (|83]) 
for Z„ at the end of inflation, we can solve (|85|) for k <C H: 



a k 



± lk/ H f£H 

Y / 7r(i — cot vir) 



k 

2H 

k 

2H 



2i 



2i 



k 

2H 

k v 

2H ) 



0[(k/2H)i 



(86) 



(87) 



Note that the above Bogoliubov coefficients (j^) satisfy the relation \a k \ 2 — \/3 k \ 2 = 1. 

We will consider the case when v = (H 2 — 4A-f§) 1 / 2 /2iJ > 0, i.e. v is real and positive, which 
corresponds to Mz < H/2. In the limit Mz — ■> H/2 (i.e. v — ► 0) conformal invariance is recovered 
and there is no particle production.^ Thus, from ([86]), to lowest order in k/H, we obtain, 



1/3* 



16tt 



2iJ 



2i/+l 



1 

V 

2 



In order to compute the spectrum of created particles, we need to evaluate, 



l^l 2 -|^)| 2 = ^(Kl 2 -i) 

l^l 2 = ^l^l 2 (89) 

where we have subtracted the vacuum contribution. In view of ( |S8|) we see that the spectrum 
of | Z | 2 on super horizon scales has a slope by k~ 2u ~ 1 enhanced in comparison to the vacuum 

9 This is because, for Mz > H/2 the quantum fluctuations of Z are damped before reaching the horizon. 
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spectrum. The enhancement is stronger for small masses Mz <C H, for which v ~ i — (Mz/H) 2 
and hence, 

where the prefactor of 2|/3fc| 2 actually appeared as (|afc| 2 + |/3fc| 2 ) — 1. 

Therefore, Z lms {l) cx I"- 1 / 2 « £-(.M z /Hf ( see also app endix |]). Thus, for M z < iT/2, the 
resulting rms spectrum for Z is almost scale invariant. Due to the (y — \) 2 factor in (pa) we 
see that the enhancement of the spectrum would be canceled if the gauge field would have been 
exactly massless, in which case v = 5. This, for instance, is the case of the massless photon. 
However, if < Mz <C H then the fact that \v — ^| ~ (Mz/H) 2 <C 1 is compensated in the above 
by considering superhorizon modes, for which H k. 

It should be noted here that in a similar way one can generate any effectively massless^] gauge 
field during inflation and obtain a similar superhorizon spectrum. A case of particular interest 
is the possibility of applying the above mechanism directly on the photon field by introducing a 
photon mass term which vanishes at the end of inflation pj| . Another realization of this scenario 
is by considering a negative coupling Hybrid Inflationary model such that the expectation value 
of the scalar field, which is coupled to the inflaton, is non-zero during inflation and becomes zero 
at the end of it. If one considers this scalar field to be coupled to the photon then the above 
mechanism is operative pjj ]. However, we feel that using the Z-field is rather more natural since, 
in this case, one does not require any additional scalars or couplings apart from the inflaton and 
the SM fields. 



4.3 The origin and nature of M z during inflation 



In computing the spectrum we have assumed that conformal invariance is broken by a mass term 
which is due to a condensate of the EW Higgs field. Here we look in more detail at the underlying 
mechanism behind the formation of such a condensate. To this purpose we need to recast (^) in 
the Hartree approximation, which we describe in Appendix |A|. We have, 



a 



Vk = 



(91) 



d 2 y k + [3A,(y 2 ), - a 2 M 2 EW - g 2 a 2 ]T ((W j ) 2 ) v - g 2 z a 2 (Z 2 ), 

j=+- 

where g 2 a 2 {-) v denotes the Hartree back reaction. Assuming that {^^) v = (y 2 ) v /a 2 , {(W^) 2 ) v 
and {Z 2 ) v are slowly varying functions in time, during inflation ( |9l| ) can be recast in the form of 
(f[8|) with, 



2 _ 9 



1 

Jr 2 



EW 



9 2 z{Z 2 



(92) 



3=+,- 



where we used d 2 a/a = 2/r 2 . As is shown in A. 2.1 , the Hartree backreaction for the gauge 
fields gives (Z 2 ) v ~ {{W^) 2 ) v < H 2 . Thus, the backreaction of the gauge fields is less than 
(2g 2 + g%)H 2 < H 2 . Also M EW <C H 2 . Therefore, the dominant contribution in ([92]) must come 
from the 3A*(^\[A term. Now assuming that v* changes slowly (adiabatically) the solution of 
© is given by @, 



,<j) 



(kr) = -V '-nrV H ( j} '(-fcr) 



C? = 1,2), 



where we used again the following normalization [cf. (|si"l)1 

W[yi 1 J(kT),yW(kT)] = Vi 



(93) 



(94) 



Assuming that Re[^*] > 1/2 then, in analogy with 



(|93|) implies, 



-"i.e. with mass < H . 
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7T 2fc 



2 

fer 



2i/«-l 



which, when subtracting the vacuum in the Hartree approximation and considering 
results in an almost scale invariant spectrum because (See A. 2. 2 ), 



ff 1 l r * 

i^ ln| 7 



(95) 
3/2, 



(96) 



where r, denotes the beginning of inflatiorj^J. The relevant quantity for the backreaction is [cf. 
(p2|)] 3A*(\l/^\l/)/7J 2 oc — ln|r|, which is a very weak function of time and thus, the adiabaticity 
assumption is justified. Moreover, for A* < 1, we see fro m (|92j ) that indeed ~ 3/2 as assumed. 
In view of (pq), the above considerations suggest that,F 



Mi 



9% 



(97) 



4.4 The generated hypermagnetic field 

Because reheating introduces large temperature corrections to the scalar potential ( |50| ) the EW- 
symmetry gets restored and the hypercharge field becomes massless. From ( |55| ) and dslf ) we 
obtain, 



(98) 



However, in contrast to the Z-boson case, the superhorizon spectrum of the photon is not 
almost scale invariant but, in fact, scales as A rms (£) <x i , because conformal invariance is re- 
tained for the photon during inflation. Thus, the contribution of the photon to the superhorizon 
spectrum of the hypercharge field is negligible compared to the one of the Z-boson. Therefore, 
the hypercharge superhorizon spectrum at the end of inflation is also almost scale invariant, 



\Y„ 



sm Ow 



6w\Z n 



(99) 



The hypercharge field during the radiation era is massless and Abelian and, therefore, it sat- 
isfies the equivalent of Maxwell's equations. The associated hypermagnetic field is defined as 
B Y = V X Y. As estimated in the Appendix the rms value over a given superhorizon scale I 
of the hypermagnetic field due to the superhorizon spectrum of the Z-boson field is, 



18(sin< 



dk 



\Z{k)\' 
V 



Inserting (|90|) into the above and after some algebra (see appendix |B|) we obtain, 



-B rm s('end) 



3\/2sin(9 w /Mz\2 



V H e J 



^end-^end 



(100) 



(101) 



where Hg is the Hubble parameter at the time the relevant scale exits the Horizon during inflation, 
kend = fc(^cnd) = 27T / £ is the physical momentum at the end of inflation, because in ( |75| ) the scale 
factor a has been normalized at that time andp| flend = H(t en d)- 



11 Remember th at, during inflation a oc 1/\t\ 

12 Considering we see that, M z ~ gz(H /2it)\J~KN where AN = JVtot - N and we used a = e N , with N being 
the number of e-foldings remaining until the end of inflation and iVtot being the total number of e-foldings. Thus (J97I) 
may be understood roughly as follows: Every e-folding the gravitationally generated fluctuation over the horizon 
volume of the EW-Higgs field is of the order of the Gibbons-Hawking temperature difik ~ H/2n. The quantity 
(\]/T\]>) represents an accumulated "memory" of these fluctuations which "pile up" while inflation continues. Thus, 
due to the stochastic nature of the fluctuations they perform a random walk, which, after AN steps is responsible 
for the V AN increase. 

13 For pure de-Sitter inflation we have Hg = ff en j. 
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Fi gure 1: Comoving magnetic field spectra and relevant seed field bounds, plotted against the comoving scale £ c . 
The steepest line (dash-dot-dot-dot) depicts the (unsubtracted) vacuum spectrum B V a C (£) oc £~ 2 . At the comoving 
scale £ c = 10 kpc one has B(;°™(10 kpc) ~ 10~ 54 Gauss. The line of intermediate steepness (dash-dot) corresponds 
to the thermal spectrum B til (£) cc £~ 3 / 2 , which meets the vacuum spectrum at the comoving scale tjj corre- 
sponding to the Horizon at the end of inflation. The solid line is the primo rdia l magnetic field, resulting from 
inflation through our mechanism, Brms(^) oc t" ~ 3 / 2 K, £~ 1 , as calculated by (105). On the scale i c = 10 kpc we 
find _B^?(10 kpc) ~ 10~ 34 Gauss, i.e. about twenty orders of magnitude stronger than the value corresponding to 
the vacuum spectrum. This is to be compared with the dynamo bounds rescaled by a factor 10 -4 , corresponding to 
the collapse-amplification enhancement factor and the factor due to the scaling between galaxy formation and the 
present (see main text). These bounds read S see d > 10 — 27 Gauss for a Universe with critical matter density and 
B secd > 10 — 34 Gauss for a Universe dominated by dark energy (fiat, low-density Universe). We also show (dotted 
line) the spectrum enhanced by helical turbulence (at £ c = 10 kpc an enhancement of about 20 is obtained). Finally, 
we note that we expect our primordial magnetic field spectrum to depart from the £~ 1 scaling law near the scale 
ill an d approach approximately the thermal spectrum since v — * towards the end of inflation. In that way the 
spectrum of our magnetic field joins the vacuum spectrum at in- 



4.5 Evolution of the magnetic field 

The creation of a thermal bath of SM particles at reheating freezes the hypermagnetic field into the 
reheated plasma in a way analogous to regular magnetic fields. At this point we should mention 
that, apart from the hypercharge component the VF-boson fields also have an almost scale invariant 
spectrum over supcrhorizon scales at the end of inflation, generated gravitationally because their 
conformal invariance is broken during inflation in the same manner as in the case of the Z-boson. 
After reheating, however, because the W gauge fields are non-Abelian they are screened and decay 
due to the development of a magnetic mass p4| , Mb ~ 0.28 g 2 T. Because of this we will focus on 
the evolution of the hypercharge magnetic field configuration. 

The physical momentum scales as k oc a -1 . Thus, by scaling k cn d to the present day we find, 

fccnd=^(^) (102) 

where t c is the corresponding scale of the mode in question today, T ren is the temperature resulting 
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from prompt reheating, T CMB is the temperature of the CMBR at present and we have used that 
a oc 7 1-1 at all times. 

If we assume that the field remained frozen until today then, because flux conservation requires 
B oc a~ 2 , we find that the magnitude of the magnetic field today would have been, 

B^ = cos 6 W BL(icod) (^) 2 (103) 

rch 

where cos 6^ is introduced due to the projection of the hypercharge onto the photon at the 
electrowea k tra nsitio n ac cording to (|56|). 



Using ( |l0l| ) and (|102j ) the above gives, 

B^Z = J(J^) 1 / 4 sin(2^)(^)(^)(^) 2 (104) 
4 407r z t c mp Hg 

where, for prompt reheating, V cno \ — (it 2 /30)g*T^ rh , with g* — 106.75 being the number of rela- 



tivistic degrees of freedom at the end of inflation.^ Putting the numbers in (104) one finds, 

bs - 5 x 10-33 (i|E£) 2 V J1 Gauss (105 ) 



where T CMB = 2.37 x 10~ 13 GcV and sia.(20w) ~ 0.84. The above corresponds to the spectrum 
of the primordial magnetic field as it would have been today were there no galactic collapse and 
subsequent dynamo amplification. Typically, such a field is referred to as "comoving" . The comov- 
ing spectrum of our primordial field, as given by (105) is shown in Fig. ^, where the substantial 



amplification compared with the vacuum spectrum is apparent. 

However, the actual physical field, being frozen into the plasma, will be affected by the grav- 
itational collapse during structure formation. Since we are interested in seeding the dynamo, we 



may use (104) to estimate the seed field at the time of galaxy formation. Scaling back the co- 
moving field to galaxy formation provides an amplification factor of (1 + z g f) 2 , where z g { is the 
redshift that corresponds to galaxy formation. This is due to the expansion of the Universe (viewed 
backwards). Moreover, the collapse of matter into galaxies brings about a further amplification 
of magnetic fields by a factor given by the fraction of the galactic matter density to the matter 
density of the Universe at galaxy formation, (pgai/pgf) 2 ^ 3 , where p g { = po(l + z g {) 3 with po being 
the matter density of the Universe at present. The above result in an overall amplification factor 
of about ~ 10 . Taking this into account and considering that the scale of the largest turbulent 
eddy corresponds to the comoving scale of about l c ~ 10 kpc before the gravitational collapse of 
the protogalaxy, we find, 



B seed ~ 10 _30 Gauss (106) 

where we assumed GUT-scale inflation Vend ~ 10 16 GeV and we have used that, for total number 
of inflationary e-foldings ~ 100 we have, (^f) 2 ~ 0.01. The above seed field is sufficient to trigger 
the galactic dynamo in the case of a spatially-flat, dark-energy (e.g. a cosmological constant or 
quintessence) dominated Universe jL0[ . 

A supplementary increase in field strength is obtained if we assume that the magnetic field 
does not freeze into the plasma upon creation, but rather that its correlation length grows quicker 
than the scale factor, as is the case for helical turbulence p5| . Such a causal mechanism can 
only operate on a given comoving scale after this scale has reentered the horizon. One can show 
that the growth of correlations due to turbulent evolution leads to an additional amplification of 
about (^cq/^e) 2 ^ 1- ^ 3 , where £ eq ~ 50 Mpc denotes the equal matter-radiation horizon today. For 
£ c = 10 kpc and v ~ 1/2 the amplification factor is about 20. 

In the above we have assumed prompt reheating. However, it is possible that reheating may 
be rather inefficient and, therefore, a long-lasting process, especially in the case when the inflaton 

14 Note that g* can be somewhat larger in simple extensions of the SM. For example in the Minimal Supersym- 
metric SM, g* = 229. Thus, g* = 106.75 may be viewed as a lower bound to the actual value of g«. 
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A 



Figure 2: The formation of a magnetic field from the generation of the Z-boson field has to pass through the 
intermediate hypercharge stage similarly to light when crossing a set of orthogonal polarizers, which needs an 
intermediate third polarizer at some angle 8 in order to go through. 



decays prominently into bosons |4q| . Indeed, the reheating temperature is typically given by, 
T Ic h < \/r in p , where T is the perturbative decay rate of the inflaton |46[||ft7|| . If the inflaton 
decays into bosons then T decreases with time and it is possible that the decay of the inflaton 
particles (inflatons) ceases before being completed leaving a small fraction of the inflatons as dark 
matter [jig fllTjl . In such case T ro h may be rather low and can satisfy the gravitino constraint 
(Troh < 10 9 GeV) even for GUT-scale inflation. Of course, an initial stage of preheating may 
increase substantially the overall efficiency of reheating ]4c|. However, in SUSY-HI, since A ~ h, 
preheating is chaotic and likely to be also rather inefficient [}I9| . Finally, one can think of variations 
similar to "smooth Hybrid Inflation" of [ |35|] , which may achieve a rather low reheating temperature. 
Still, as shown in (101), the amplitude of our hypermagnetic field is, in fact, decided by the 
inflationary energy scale and not by the reheating temperature. Indeed, it is easy to show that, even 
in the case of long- lasting, inefficient, reheating the strength of our seed field at galaxy formation 
remains essentially unmodified, so that ( |105| ) is still applicable. Therefore, our mechanism may 
manage to satisfy the gravitino overproduction constraint while still able to provide a sufficiently 
strong seed field. Thus, our SUSY-HI inflationary model can be incorporated in a Supergravity 
framework without problem. 

At first sight, it may seem unlikely that a magnetic field is obtained from the gravitational 
production of the Z-boson field, which is orthogonal to the photon field. In fact, the situation is 
analogous to that of light polarizers (see Fig. ^). With two orthogonal light polarizers no light 
passes. When a third polarizer is inserted at an angle 9 (with respect to the second polarizer) 
however, some of the photons do pass. The photon amplitude is reduced by ^sin(2#), just as in 
the case of the Z-field. The main advantage of the amplification mechanism presented here is its 
naturalness. Indeed, no fields are required except those of the standard model and the inflaton. 
Moreover, the mechanism is independent of the particular model of inflation considered, but can 
be thought of as a generic consequence of inflationary theory itselfpl 

An exciting possibility is that the above magnetic field can be amplified even further during 

15 Note that most of the successful inflationary models (e.g. chaotic, natural, hybrid) have GUT-scale V en d- 
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preheating in an explosive resonant way. In what follows we explore this possibility using the 
particular model of Flipped SU(5) SUSY-HI. 



5 Preheating 

In addition to the growth of the superhorizon spectrum due to the breaking of the Z-conformal 
invariance in the inflationary period, the generated hypermagnetic field may be amplified more 
during preheating due to the backreaction of the hypercharge source current. From (ftp| ) it is 
evident that the growth of the amplitudes of the GUT-massive gauge fields and Y" is a 
prerequisite for this source current to be substantial. In Hybrid Inflation this is achieved through 
their mass term given by which resonantly pumps energy from the GUT-Higgs field. The 

latter is resonantly produced at the end of inflation due to its coupling with the oscillating inflaton. 



5.1 Parametric resonance in Hybrid Inflation 
5.1.1 The mode equations of the scalar fields 

We shall consider production of the supermassive vector bosons by parametric resonance caused by 
oscillations of s and <j>. Since s is a singlet and does not couple to the vector fields, the production 
of gauge fields is due to the GUT-Higgs field. A significant production of modes of any kind occurs 
only after <fi and s start oscillating around the bottom of the potential. Oscillations of cj> are not 
harmonic, but to a good approximation they are periodic, with an almost constant period. This is 
so since the Hubble constant is small, H = ^/{p- m f/2>)/mp ~ XM 2 /mp ~ 10 13 GeV, while a typical 
oscillation frequency is ui ~ y/XM, where we have used ([!]) to estimate p- m { = V(0, s c ). 

Since, typically, the resonant decay of a field lasts a few hundreds of oscillations, the decay 
time of the inflaton is < H . On these grounds one can neglect the expansion of the Universe 
during preheating and regard the scale factor to be approximately constant a ~ 1 in agreement 
with the normalization (f75j). Therefore, the field equations ( p!5| ) and ( |4S| ) become, 



= 



d 2 + \(<p 2 - M 2 ) + ha 2 ip 







(107) 
(108) 



where in (108) we have ignored the backreaction of the supermassive GUT-bosons because at the 
onset of preheating their amplitude is negligible. In terms of the field (p — ip — M the above can 
be rewritten as, 



d 2 + A (Cp 2 + 3<pM + 2M 2 ) + ha 2 



(p + ha 2 M 



d 2 + h(p 2 + 2(pM + M 2 ) a = 



(109) 



where we have ignored the term since the slope of the V a potential is negligible compared to 
the interaction term. The solutions to these equations are damped periodic oscillatory functions^] 
which in general have the following form, 



cp = tpoP v (v<pT) 

a = ff P a (u„T), (110) 

16 We do not consider here the effects of chaotic resonant behavior which may occur in hybrid inflationary 
models. 
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where P ip (x + 2tt) = P v (x), and P a (x + 2ir) = P a (x), maxjlP^I} = 1 and max-dP^I} = 1, and (po 
and (To are the amplitudes of oscillations, which are in principle time-dependent since ip and a 
eventually decay and redshift. 

The corresponding mode equations can be obtained with Fourier transform. The result is, 



d 2 + k 2 + A (3(p 2 + 6<pM + 2M 2 ) + ha 2 <p k + 2ha(<p + M)a k = 
d 2 + k 2 + h{(p 2 + 2q>M + M 2 )} a k + 2ha{(p + M)ip k = 



(111) 
(112) 



By making use of (<p 2 ) = (fil/2 and (a 2 ) = (7q/2, we can easily estimate the effective masses of 
the modes to be, 



nr = 2A ( M 2 + ^ J + i ha 2 ~ 4AM 2 



in' 



h M 2 



- 2 



-hlVP 



so that the dispersion relations read, 



(113) 



k 2 + 2\ { M 2 + -^>l 



-ha 2 
2 



k 2 + 4AAf 2 



-KM 



(114) 



where ip\ ~ M and a 2 ~ s 2 — (X/h)M 2 . Note that the effective mass of the ip modes is a factor 
of \/2 larger than the Higgs mass Mh — V2XM. 



We need to investigate whether the mode equations ( 111 ) and ( 112 ) result in a broad or narrow 
resonance, and also whether the leading (broad) resonance is active. To do that we have to consider 
the quality factors q for the mode equations and then compare them to the effective masses (113). 



5.1.2 The characteristics of the resonance 

The generic form for a resonance equation is 



d 2 4> k 
d 2 zu 



A k + S ^2q l P i (-cu) 



V^fc = 



where 



Ai 



m^g + k 2 



(115) 



and qt are the quality factors with oscillating functions Pi{w) ((fi) ro = 0, max^jPi} = 1) whose 
period is either 7r or 27r, and w = lo v t is the rescaled time variable. In the Mathieu case, 
Pi = cos(2nj). Some of the references where one can find instability charts, which are plots of 
the Floquet exponent /i in terms of q an d A, are listed in [^o| jl?) . 
Consider first the quality factors for (111), 



3A^| 
4< 



3 
16 



lip 2 



6XMp 
2uj 2 



3 
4 



2hMcr 



(116) 



Since uj^ ~ 4AM 2 , we conclude that all the quality factors of (111) are generically of the order 
unity. The quality factors for (112) are on the other hand 
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h<pl 1 h 2hM(p 1 h , 2hMa 1 h 

9ct1 = 4^~16A 9CT2 = ^r^4A ^ 1 = ^T^4Va (U7) 



so that they are typically smaller than the quality factors (116) of l \l 1 1| ) .P| Further we can compute 



ml . 4 



A v (k = 0) = -f = l~-q v2 



u% 3 

A a (k = 0) = <^ = l^ qa2 (118) 

which then determines onto which resonance the field can decay. Note that since the position 
of the first resonance is at A ~ 1 for q < 1, we conclude that the <p field tends to decay into the 
first resonance such that the infrared super horizon modes get populated (see figure in |4(j]) with 
fi ~ 0.3 or smaller. A similar statement applies to the a field. 

The strength of the broad resonances may be as large as /i ~ 0.2 along the A = 2q line, and 
becomes rapidly stronger at smaller values of A, i.e. A <2q (cf. Eq. (7) in ]5Cfl). Even if the field 
can decay only into the second available resonance, the strength is expected to be of the order 
/12 ~ 0.04, just like it is in the case of the chaotic inflation with a quartic scalar interaction term, 
which decays into the second resonance. 

Regarding the a field, depending on the choice of the couplings, q a i can be larger or smaller 
then unity. Generically, however, we expect that A ~ h, and hence q a 2 is typically of the order 
(but smaller then) unity. In SUSY-HI we have A = 2h and, therefore, q a 2 — g < 1. 

In the above, by taking a ~ 1 we have ignored Hubble damping. It is known that the latter, in 
fact, results in stochastic resonant amplification, in the sense that the action of the damping term 
may be positive or negative on occasion, depending on whether the amplitude of the oscillating 
fields is increasing or decreasing. However, as has been shown in p7[||i^| , the net effect of such a 
term is, in fact, in favor of resonant production. 



This analysis has not covered the mixing terms in ( |lll| ) and (112). Numerical simulations 
done by the authors of ]50| indicate that these terms typically cannot prevent resonant decay, and 
in some cases can enhance it. 

5.2 Resonant production of massive vector bosons 

The field equations for the supermassive GUT-bosons are ( |23| ) and (p4]), while j2^ ) is the one of 
for the U-boson. The resonant production of the massive bosons is due to the fact that their 
equations of motion are of harmonic form with frequencies given by their oscillating mass term. 

5.2.1 The parametric resonance equations 

At the onset of preheating we may ignore the backreaction source current because the initial 
amplitude of the massive bosons is negligible. Later on the source current contributes additionally 
to the production of the massive bosons. 



Let us concentrate on (|23|) first. Proceeding in a similar manner as in 4.1 we obtain for the 
longitudinal and transverse modes, 

^ + k- + Mhr)+ 2 %\ M ^ d T ]?r\\ = (119) 

2 ^ v a± 

where 



Mj CY )X aM = (120) 



17 Remember that in SUSY-HI A = 2h. 
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/cPk 



(121) 



Recalling that ip = M + lpqP v {uo v t) and using ( |2l| ) we can find the effective mass A4 XY of X 
during preheating to be, 



M 2 XY = -g 2 (M 2 + + 2M l p P lp ) 



M 2 XY = -g 2 M 2 



2 V M ; 



-g 2 M 2 



(122) 



wher e we used that (P 2 ) T = 1/2 and (P V ) T = 0. Rescaling the time as uj^t —> ra, the equations 



(119) and (12C) can be recast as, 



k 2 +M XY \ 2(d VJ \nM XY )k 2 



k 2 + M 2 



XY 



M 



XY 



X a " = 



X 



= 



(123) 
(124) 



Now, it easy to show that, in the above, we have, 



k 2 + M 2 XY 



k 2 +M\ Y . 5 2 M 2 r l t^o 2 ~ 



2^ 2 



TM' 



2^P V 



(125) 



where P v — 2P 2 — 1 is a periodic function of w with frequency ui v = 2uj v and max{|P v |}=l. Each 
of the above equations contains two resonant channels with the g-factors, 



Qi 



g 2 M 2 tpo g 2 



2lo 2 

***** <p 



M 8A 



g 2 M 2 yp 2 ^ g 2 _ 1 
%ujI { M> ~ 32 A ~ 4 91 



M 2 XY 



16A 



-qi ~ 6q 2 



(126) 



This implies that, provided q\p > 1/2, the q\ -resonance is efficient in amplifying gauge fields, 
while the qi -resonance is inefficient. In fact when Ak=o < 2q {Ak=o > 2q] the growth rate of the 
field grows {decreases} with increasing q |5Cj], implying that for Ak=o — §9i the resonance becomes 
stronger as q\ increases. For example, one enters the first broad resonance at about q\ ~ 1/2 and 
exits it at about q\ ~ 3. 

When compared with the transverse equation (|l24| ), the longitudinal equation ( |l23| ) contains 
an additional 'damping' term. The sign of the 'damping' coefficient d^hiM XY = d^,P v /P v may 
be either positive or negative, resulting in damping or growth, respectively. Since the positive and 
negative parts are symmetrically and evenly distributed, we expect that the growth wins over the 
damping, just as it is in the case of stochastic resonance p8[ . As a consequence, the longitudinal 
modes grow faster than the transverse modes. Large longitudinal amplitudes can be mediated to 
the transverse modes through scatterings. We do not consider here possible physical consequences 
of these processes. 

Similar results are obtained for the Y° and the bosons, substituting, in the latter case 
Mxy — ► My or equivalently g — > 2^/ (6/5) gj cos0. 
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5.2.2 The amplitudes of the massive vector bosons 

In order for the backreaction to stop resonant production of the massive gauge fields, the induced 
(Hartree) mass must be larger than about the mass of the Higgs boson, since this is the typical 
resonance scale. This can be argued as follows. The shift in A required to switch off the leading 
resonance is of the order 5A ~ 0.5. This cannot be mediated by the (^-field modes, since the 
maximum possible expectation value is ((<5</9) 2 ) max ~ M 2 , resulting in S A = 5M 2 XY /oA ~ g 2 /16A, 
which is typically not sufficient to shut down the resonance. Moreover, the inflaton does not help, 
since it does not couple to the vector bosons. 



There are however cubic terms in the source current, which have been neglected in (124) and 
which may shut down the resonance. The cubic source current terms can be obtained from (|2F 
Similarily to (|32|), for the supermassive GUT-bosons we have, 



(Wy C = g^g^^mM^phi^jc, - (p <- a)} (127) 
where £ = 4, 5 and a = 1, 2, 3. In view of (]13]), (127) can be written schematically as ^| 



Jxy ~ 9 {XXX + WWX + WVX + VVX), (128) 

where X denotes the supermassive GUT-bosons X^ and Y" W denotes the electroweak gauge 
bosons Wp, Yp, which are gravitationally produced in inflation as shown in Sec. ^, and V is the 
massive U-boson. Note that the contribution from the G-bosons of SU(3) C is negligible because 
they are neither produced gravitationally in inflation, nor are they generated during preheating, 
since they do not couple to any Higgs field. Thus, a typical contribution of the backreaction 
current to the mass of the X® Y° fields, coming from their resonant production, is of the form 

6M 2 XY ~ g 2 ((X 2 ) + (V 2 )) (129) 

5 



Similarly, for we have from ( |28| ) that Jy = ~ \ Ti cos ©^y so that the cubic term contri- 



bution is given by (|34| ) to be, 



Jv h[c ~ I g 2 XXV (130) 
Hence, the contribution of the backreaction current to the mass of the is, 

6My~lg 2 (X 2 ) (131) 

D 

In order to get the first resonance to shut down, one requires SM 2 ~ Mfj — 2AM 2 . Hence, the 
energy stored in the X and V fields is then of the order, 

Pxv - 9 2 ((X 2 ) 2 + (V 2 ) 2 ) ~ ^M 4 ~ A pinf (132) 

9 9 

Since A ~ 1 we conclude that the resonance shuts down only when most of the energy is in the 
massive gauge field modes. In other words, the resonance is efficient until a significant fraction 
of the inflaton's energy decays. Then, as the oscillating amplitude of ip decreases, one enters the 
narrow resonance and the decay slows down. This is just like the case of two scalar fields with 
q ~ 1, and also the one field case, where also q ~ 1. One is then left with the condensates of the 
gauge fields whose amplitude can be estimated to be, 

\X?\, Kl K\ - 9* 1/4 — - 0.1M H (133) 

g 

where g* < 100 is the effective number of degrees of freedom produced by the resonance. 

18 Here we have used that — (— — h c' ) = 1. 
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5.3 The hypercharge field at preheating 

In a similar manner as for the supermassive GUT-bosons one finds from ( p9| ) that the hypercharge 
transverse mode equation is, 



(d 2 T + k 2 ^y^ = y^sinej-^ 



(134) 



where 



/d k f d k 

jy-^J^(k,T)exp{ik • x) and Y v (x, t) =J ^— ^ y" (fc, r) exp(ifc • x) (135) 

and we have used the current conservation equation, d T {Jfy)T = ik • 3y . 

Using t he H artree approximation, as shown in the appendix |A.l| , and also in view of (^) we 
can recast (134) as, 



{d 2 + k 2 + 6f Y [l(\x\\ - (\x T { 2 ) v ] } y 1 - = -efycot e [l<|X| 2 ) v - (\x T 



(136) 



where g 2 , = Iffy- and X^ stands for all the 6 the supermassive GUT-bosons X", Y" with, 



1 f k 2 dk 

(\Xt\ 2 ) v = v J ^-X T (k,r)X T (-k,r) 
and also [cf. <M) and @], 



1 

r 
i 

r 



k 2 dk 



2vr 2 



2tt 



|A-(fc,r)| 2 
a [^r(fc,r)] 2 



(137) 



(l*r| 



1 

r 



fc 2 [a r Af"(r)] 2 
fc 2 + M| F 



(138) 



We shall now investigate whether the hypercharge field ( |136| ) may undergo resonant amplifica- 
tion at preheating. Since the supermassive GUT-bosons are mainly produced at preheating, the 
dominant contribution to X 2 ms comes from the resonant modes that initially oscillate in phase |30) , 
which in principle may drive resonant amplification of the hypercharge. We now discuss the nec- 
essary conditions required for this indirect resonance to be operative. 

At early times in preheating the Hartree terms X 2 ms are small, so that the hypercharge bosons 
grow through (inefficient) narrow parametric resonance |5l|| . When the oscillatory component of 
X 2 ms becomes of the order ui 2 however, the hypercharge resonance can become broad Since 
the hypercharge field has an infrared spectrum that is already amplified in inflation, while the 
U-boson spectrum is that of the Minkowski vacuum, when compared with the X 2 V-term, the 
X 2 Y-tevm is expected to dominate the resonant growth on superhorizon scales. Therefore, we 
may write, 



J$(t) -6V2gg Y X 2 ms y x (T) 



so that (136) simplifies to 



2-M-L 



d 2 y 



dr 2 



(k 2 + 6g 2 Y X 2 ms ) y x ~0 



(139) 



(140) 



Now, X 2 ms can be conveniently split into the slowly growing and oscillatory contributions as 



follows: 



X 2 ms = X 2 (r)+X 2 sc [l + P x (r) 



(141) 
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where Px (t) is a periodic function with period ~ tt/uj^ such that (Px(t)) t — and (P x (t)) t = 1/2, 
and 

X 2 (r) = A 2 ac (e 2 ^ - 1) , A 2 sc (r) = A 2 ac (e 2 ^—- - l) (142) 

where [ix denotes the relevant resonant growth rate of the supermassive GUT-boson fields X, and 
Ho is the growth rate of the constant contribution to the A-field mass induced by the inelastic 
scatterings. Also, we have subtracted the vacuum contribution X 2 ac ~ fc 2 os , with fc 2 0S being the 
typical resonant momentum. Provided fig < MX; or equivalently Xq < X osc , the hypercharge grows 
resonantly. The amount of available energy p ln f ~ AM 4 /4 results in the following upper bound, 



~M 2 XY X? ms < ~M 2 H M 2 (143) 



or equivalently [cf. (|122|)] 



■Y.i,,- <I^ 2 M 2 (141) 



where .M^y ~ 3g 2 M 2 /4 and Mj^ = 2AA/ 2 . This is in agreement with (133) and implies the 
following bound for the quality factor of the hypercharge resonance, 



2 



= 6 9yX < __ sin , @ < 4 (145) 
2w£ 12 

where sin0 = gy/g. Since significant resonant amplification of the hypercharge on superhorizon 
scales may occur without tuning only when qy > 1/2, we conclude that the inflationary spectrum 
of the hypercharge gauge bosons is hardly amplified by parametric resonance. 



6 Conclusions 

We have presented a mechanism of primordial magnetic field generation based on the breaking of 
conformal invariance of the Z-boson field during inflation. The mechanism is generic and inde- 
pendent of inflationary model as long as the reheating temperature is higher than the electroweak 
scale. This is because our mechanism requires a phase of electroweak unification to "channel" the 
generated superhorizon Z-boson spectrum into the photon, through the hypercharge field. 

The conformal invariance of Z is naturally broken due to its standard-model coupling with 
the electroweak Higgs-field. The latter, during inflation, develops a condensate comparable to 
the one of the inflaton field itself resulting in a non-zero, but small, mass for the Z-boson field. 
Conformal invariance breakdown results in gravitational production of the Z- field on superhorizon 
scales. We have computed the relevant superhorizon spectrum and found that it is almost scale 
invariant. At the end of inflation reheating restores the electroweak symmetry and the Z-spectrum 
is converted into a hypercharge spectrum (the photon's contribution being negligible) , which, due 
to the stochastic nature of the original Z-fluctuations, gives rise to a superhorizon hypermagnetic 
field that freezes into the reheated plasma. We have calculated the spectrum of the rms value 
of the hypermagnetic field and found that it is of the form B Tms oc \ jl as shown in Fig. |l} In 
a similar way non-Abelian VF-fields are also generated but their associated magnetic fields are 
screened due to the existence of a magnetic mass of non-perturbative nature, which stems from 
their self-interaction. During the radiation era the hypermagnetic field evolves satisfying flux 
conservation. At the electroweak phase transition it is transformed into a regular magnetic field. 
When scaled until galaxy formation, this magnetic field is found to be sufficient to trigger the 
galactic dynamo and explain the observed galactic magnetic fields in the case of a spatially- flat, 
dark-energy dominated Universe with GUT-scale inflation. 

The beauty of our mechanism, apart from being independent of the inflationary model, lies in 
that it does not, unlike most proposed mechanisms, require the explicit introduction of conformal 
invariance breaking terms in the Lagrangian or any exotic fields other than the ones of the standard 
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model and the inflaton field. Moreover, we do not have to specify a particular GUT group, or 
involve grand unification in any way. 

An intriguing possibility was that our magnetic field could be further amplified during preheat- 
ing. In order to study this we considered a (supersymmetric) hybrid inflationary model and grand 
unification under Flipped SU(5), both well motivated. However, we have found that preheating 
amplification is probably negligible since the hypercharge field is amplified only via indirect narrow 
resonance which turns out to be rather inefficient. Still, a better (possibly numerical) treatment 
of the magnetic field at preheating is probably necessary to provide an adequate understanding of 
its behavior. Such treatment should, in principle, address issues such as the conductivity of the 
newly created plasma or the fate of the hyperelectric field. 

In particular, the behavior of the conductivity during preheating may inhibit our hypermag- 
netic field. However, an analytic treatment of this issue is rather complex and beyond the scope of 
the present article and numerical studies are inconclusive. It seems that the behavior of the con- 
ductivity (spatial distribution and growth rate) is highly non-trivial, non-perturbative and model 
dependent [p2f . It has been shown however, that the appearance of conductivity in preheating 
does not necessarily inhibit the growth or the existence of (hyper)magnetic fields but may allow 
substantial amplification ]2^] |^3| . This is certainly an issue that deserves further investigation. 

In summary, we have presented a natural mechanism for magnetogenesis by inflation which 
may be an explanation for galactic magnetic fields. 
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A The Hartree approximation 

In this appendix we discuss the use of the Hartree approximation in the computation of scalar 
and gauge fields. This approximation is applicable in this case mainly due to the classicallity of 
the superhorizon modes and also since the typical value of the quality factor q of the parametric 
resonance is about one or smaller, essentially until the field grows to its maximum value -p*] 

The Hartree approximation takes account of elastic scatterings only, which do not change the 
momenta of incoming particles, and it models the dynamics well when 2-to-2 scatterings are the 
only ones that are relevant. This is indeed so at the early stages of preheating, as is known 
from exact classical simulations (e.g. of the scalar theory). When thermalization begins, inelastic 
scatterings, which include 2-to-4, etc., become important. They combine the infrared modes into 
more energetic ones, as it is required by the thermalization process, because particles produced 
by parametric resonance are more infrared in comparison to a thermal distribution. The Hartree 
approximation breaks down after a significant fraction of the inflaton's energy has decayed. 

The Hartree approximation consists essentially in replacing a product of two fields (in an 
equation of motion) with the spatial average as follows, 

A(x, t) B(x, t) -> (A(x, t) B(x, t)) v = AB + (SA(x, t) SB(x, t)) v (146) 

where A and B denote the zero modes defined as A = (A(x,i)) v , B = (B(x,t)) v and also 
6A(x, t) = A(x, t) — A, 5B(x, t) = B(x, t) — B, with the volume averages being defined as, 



1 I" 1 f' ii 3 k 

(SA(x,t)SB(x,t)) v = - J d 3 x6A(x,t)5B(x,t) = - J —A(k,t) B(-k,t) 

A = — fd 3 xA(x,t) = y A(k = 0,t) (147) 



19 The validity of the Hartree approximation can be checked by performing a consistent one-loop calculate 
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In these definitions the volume V is defined by the following discretized version of the momen- 
tum integral 



f d 3 k 1 v-^ 

iw^? 11481 

The prime in the momentum integral indicates that the zero mode should be taken out of the 
integral. The zero modes require special care, since they may become macroscopic, leading to a 
condensate.^ 

A.l Hypercharge source current 

Let us employ the Hartree approximation in computing the source current of the hypercharge 
which is shown explicitly in ( |38| ) . In the spirit of our approximation, we shall neglect the terms in 
which the zero modes figure incoherently in quadratic or cubic combinations. This means that at 
this point we shall take account only of those nonlinear contributions from the zero modes which 
oscillate coherently in time. 

Some of the remaining terms are composed of the fields which have incoherent phases in 
different points in space, so that, when averaged over space, they vanish. Including these terms is 
strictly speaking beyond the Hartree approximation^ Therefore, to first approximation we keep 
only the terms that contribute to the genuine Hartree approximation, which are the terms that 
oscillate coherently in space and time, since they dominate the resonant production. 

Consider first the derivative terms of the source current in (j38[), 

- ^ sin 9 r Y [XdX] = - g Y Im{ + In ^D- g )]g»0g™(X?X a a + Y?T* a ) + 

(149) 

+ g^g^[X»d p x" + Y«d p Yl " (P ~ *)]} 
In CFRW with a ~ 1 the above may be written as, 



'A sin0 jJ[XdX]=-s y ,rim^ 

(150) 

Hartree averaged the terms in the above are of the form, 



{XZdiX 



y { d\X«d X a a = ' 



d 3 k 



(151) 



V 



d 3 xX^X a 



d 3 k 
V J (2tt) 



3 M 



where Xff is defined in (121). As a consequence of temporal and spatial isotropy, we have 
(X«d p xZ)v = X™dpC « »7pa that is, 



(X%d p X a a ) v = l%«d p J2 \X«\ 2 => Im((JT«flX)v) = 

a 



(152) 



20 Note that this ensemble averaging corresponds to the classical ensemble representation of the quantum state. 

21 One way of improving on the Hartree approximati on i s to treat the terms that incoherently contribute as noise. 
When averaged over time, these terms yield zero. In fea | it has been shown that including these noise terms can 
only enhance resonant production. Thus, when not considering these terms we adopt a conservative approach on 
the resonant growth. 
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where there is no summation over the a index. Since all the terms in (|150j) are of the same form 
we conclude that, in the Hartree approximation, jJ[XdX] = 0, i.e. the spatial averages of the 
quadratic derivative terms in (38) average to zero and only the cubic terms survive. 

Let us consider the cubic terms now. We start by noting that, for two different gauge fields 
and B p we have, (A p B v ) v — r\ pv A p B p . The zero modes are proportional to the polarization 
vectors, that is, 



iP=-^(fe = 0,r)c<e^^ and B*' = -^'(fc' = 0, r) oc e^ , (153) 

where p — T, L (T — 1, 2 and L = 3) denotes polarization, and ^p^i^yY = — ■ ^ = —5k,k' $ pp ■ 
As a consequence of the vectorial nature of the condensates, and their origin in quantum 
fluctuations, we conclude that they are randomly oriented, and hence, for any different gauge 
fields, 

{A tl B v ) v =r ltlu A p B<' = Q (154) 

In view of the above result the hypercharge source current in ( pl[ ) in the Hartree approximation 
can be recast as, 

-^J sine J Y = - g*£^Be((XpC)v + <X?K)v) y p~ 

- C ot®? Y £^Re((X<pC) v + (yffiy) V p - (155) 

Because the X p and Y" bosons are entirely equivalent we define, 



(X^X a ) v = (X°X*) V = (Y«Y a a ) v V a = 1, 2, 3 (no summation) (156) 

Thus, the above becomes, 



% sine = -6g^ pa Re^(X p X a ) v j Y p - 6 cot 6 g Y £ pvpa Re({X p X a ) v J V p (157) 



Therefore, in the Hartree approximation only the terms of the Y p and V p bosons contribute to 
the hypercharge source current. In view of the definition ( |39| ) it is straightforward to show that, 

(158) 



£^Rc({X p X a ) v ) = Re{X"X V ) v - r,P»(X a X a ) v 
Thus the source current becomes, 

- sinG J Y = -6g^(Re(XPX v ) v - ^ (X a X a ) v ^j (Y p + cot Q V p ) 
In (|136|) we need the spatial component of the above, for which we find, 



(159) 



£ sine(Jy)i = 6g 2 Y [(\X \ 2 ) v -~(\X\ 2 ) v \(Y i + cote V t ) 



(160) 



where \X \ 2 = X Q X , \X\ 2 = X 3 X } and we used (X 1 X 1 ) V = (X 2 X 2 ) V = (X 3 X 3 ) V = ±(|X| 2 ) V . 
Because (•)„ is a function of time only, we can immediately obtain the Fourier transformed current, 



j2 sin 6 J Y = 6 g Y 



l*r| 



-|<|X| 



(y 1 - +cotev ± ) 



(161) 



which is used in (136). 
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A. 2 Gravitational production during inflation 



Here we will briefly calculate the volume averages of scalar and gauge fields which are gravitation- 
ally produced during inflation. In particular we will focus on the iT-boson field and the electroweak 
Higgs field tp. 

A. 2.1 The case of the Z boson 

For the Z-boson field we have, 



1 f d^k r 

(Z 2 ) v = (0\Z ■ Z\0) (Q\(Z ■ Z) vac \Q) = - j— £ \Z^(kr)Z^{kT 



V 



2c4(fc,r) 



(162) 



where p denotes the three possible polarizations, w z {k,r) denotes the dispersion relation for Z 
and for the vacuum. 



V 



2u p z (k,T) 



1/2 



3 —i J uj p '(k,r)dr 



Now, using (fed) and taking ujz — k (U63) becomes, 



(163) 



(Z 2 



H 



k dk 
o 4^ 



Then, considering | — v ~ (Mz/H) 2 <lwe get, 



2 

kr 



(164) 



(Z 2 



-3H 2 



du 
o 8^ 



-(M z /H) 2 



3H 2 (M z 



8tt 2 V H 



where u = (kr) 2 . Obviously, one arrives at a similar result also for the VF-bosons. 



(165) 



A. 2. 2 The electroweak Higgs field 

Similarly for the EW-Higgs field one has, 



(y 2 



\y 2 \0) - (Oli&JO) 



d A k 
V J (2^f 



y (1) (kr)y w (kr) 



(2), 



V 



where u)^(k, r) denotes the dispersion relation for the scalar field and for the vacuum, 



Z/vac(fe,T) 



V 



_2w^,(fc, r) 



1/2 



-i J uJtp (k,r)dT 



Using then (95) and also u)^ ™ k, we obtain, 



H 



k dk 



2 

kr 



With v* 



| the above gives, 



a 2 H 2 



1 du_ 

8^ 



a 2 H 2 
4tt 2 



■In 



(166) 



(167) 



(168) 



(169) 



where the infrared cutoff corresponds to the onset of inflation and, for supcrhorizon scales, 
k\r\ ~ k/H < 1. The above directly results in (pKih.P 

Comparing the backreaction terms in view of (165) and ( p6| ) we see that the 3A*(* t *) v term 
is indeed the dominant when A* > g 4 /2w 2 ~ 10~ 3 . 

22 Here we used that k = aH = — 1/r => ln|fci/fc2| = hi l^/nl- 
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B The rms value of the hyper magnetic field 

Here we will calculate the rms value of the hypermagnetic field Bj ms (i) as a function of scale at 
the end of inflation. 

The definition of the hypermagnetic field implies, 

Bj{k- 1) = ie^ l k 3 y t = ie ijl kjy^ (170) 

where we used, e i: > l kjyf = and 

B Y (k,t) = J ' d\e-^ x Bj{x,t) B Y (x 7 t)= J ^e^ x B Y (k,t) (171) 

?<(*,*)= J d 3 xe- ik - x Yi(x,t) Yi(x,t) = j ^e a '^(M) (172) 



We now define, 



where 



{B Y ms f = (B Y (l,x)B Y (£,x)) ee ly d\Br{l,x)Bj{l,x) (173) 



5^,*) = ! / D)(x') (1741 



Thus, in view of the above we have, 



(B r Y ms ) 2 = 4 / dh f d?y j dh j exp[*(fc -y + k'. z)] B Y {k)-B Y {k') (175) 

v J J\K-y\<iJ\v-x\<t J y Z7T ) 



Now, since B y '= ikxy ± (k) we have, 

B Y {k)-B Y {k') = [k- y ± (k / )][k'-y ± (k)} - (k • tf)[y ± Qe)-y ± (tf)] (176) 
Using this and substituting K = = (fc + fe')j K = o(fe — and r*i = y — x 7 r 2 = z — x we 



can recast (175) as 



(BLs) 2 = ^J d 3 x Jd 3 n Jd 3 r 2 J f0^ex P { l {K.(r 1 +r 2 + 2x) + K.(r 1 -r 2 )}} x (177) 

: {K-[y x (K+K) + y ± (K-K)] - K-[y^(K+K) - y J -(K-n)] + (k 2 -K 2 )[y x {K+K,yy x {K-K)]} 

where K = \K\ and k — \k\. Since, J d 3 x J -^jje 2 ' lK ' x — %S(K) the above becomes, 

(BLs) 2 = 4~3 I d ^ I d ^ / tSj e ik <^e\y(k)\> (178) 

J \r- t l<t J lr 2 l<?- J \ > 

where we put k — ► k and^J y^^kyy^^k) ~ 4|3^(fc)| 2 (because the transverse component has two 
polarizations) while also using, k-y ± (k)= fc-^V-^-fc) = 0. 
It is easy to show that, 

23 Because the coefficients for Z 1 - in (fi 7 ^) do not depend on direction, it is reasonable to assume that, on average, 
the hypermagnetic field depends only on the magnitude of the momentum, so that y (fe) = y ' (fc). where fe = |fe|. 
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d\ e ±ik • 



2nr dr / d(cost9)e 



±ikr cos $ 



47T 



sin(fcf) -£cos(k£) 



(179) 



where cost? = (k • r)/kr and r = |r|. Using the above into (U7q) we obtain equation (100), 



18 (sin 9 



+ OC 



-sin(fcf) -£cos(k£) 
k 



l^(fc)| S 



fliog) 



where, from (|99|), y 1 " — — sin (9^-2^ and we took V = §7r£ 3 . Since we expect that the contribution 
of the modes with k > 1 is negligible we may consider only the supcrhorizon spectrum of Z^-^fc), 
which is the one given by (BOh . Inserting the latter into (100) and after some algebra we find, 



9 (sm6 w ) 2 H 2 fM z 



(2£H) 



H 



-t-oo 



dw 
w 



sin w cos w 



(180) 



where w = k£. Using spherical Bessel's functions the integral / in the above evaluates to, I = 1/4. 
Therefore, we obtain, 



Sv^sin^ (H\(M. 



H 



(181) 



2£(2£Hy 

In fact, one obtains, Bj ms {i) oc g u ~ 3 / 2 which, in the limit when Mz — > i?/2 (^ — > 0) corresponds to 
the thermal spectrum of massless particles, B ex l^l 2 . However, we actually have (Mz/H) 2 -C 1 
so that, the above becomes (101). 

In a similar way one can obtain the rms value of the Z-boson field over superhorizon scales. 
One finds, 



^2T[(Mz/H) 2 } ( H ]( M z] 2 {Mz/H? 



(182) 



8(2^i/) 1 

We see that the spectrum is almost scale invariant. Note, that the value of the Gamma function 
is very large for (Mz/H) 2 <C 1. Physically, this is because the scale invariance of the spectrum 
results in all the modes with k < l/£ contributing substantially to the rms value. 
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